Fractionalization in a square-lattice model with time-reversal symmetry 
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We propose a two-dimensional time-reversal invariant system of essentially non-interacting electrons on a 
square lattice that exhibits configurations with fractional charges ±e/2. These are vortex-like topological defects 
in the dimerization order parameter describing spatial modulation in the electron hopping amplitudes. Charge 
fractionalization is established by a simple counting argument, analytical calculation within the effective low- 
energy theory, and by an exact numerical diagonalization of the lattice Hamiltonian. We comment on the 
exchange statistics of fractional charges and possible realizations of the system. 



Introduction. — It is now well-known that fractional quan- 
tum numbers can arise as the collective excitations of a many- 
body system. The canonical example of such fractionaliza- 
tion is a 2D electron gas (2DEG) placed in a transverse mag- 
netic field in the fractional quantum Hall regime. At odd 
inverse filling factors, v > 1, the many -body ground state 
is described by a strongly-correlated Laughlin wave function 
and the time-reversal symmetry is broken. The excitations 
carry the fractional charge ve (Ref. 1) and exhibit fractional 
(Abelian) statistics,- The search for other systems that ex- 
hibit fractionalization is ongoing. Two important questions in 
this search are whether strong correlations or a broken time- 
reversal symmetry is necessary for fractionalization to hap- 
pen. 

Recently the answer to both questions is argued to be neg- 
ative. A group including the present authors 3 proposed a sys- 
tem with fractionally charged, anyonic excitations that can be 
described by a weakly-interacting wave function found by fill- 
ing a set of single-particle states. Moreover, Hou, Chamon 
and Mudry 4 have argued that a vortex in the Kekule modula- 
tions of the hopping amplitudes on a honeycomb lattice, like 
that of graphene, binds a fractional charge e/2 without break- 
ing the time resversal symmetry. 

In this paper we propose a system on a square lattice with 
time-reversal symmetry that exhibits fractionalization. The 
system consists of a square lattice threaded by one half of a 
magnetic flux quantum cDo = hc/e per plaquette on which 
electrons can hop to nearest-neighbor sites with no interaction. 
Time-reversal symmetry is preserved in this lattice model be- 
cause electrons cannot detect the sign of the flux with mag- 
nitude <I>o/2. In addition, we assume a dimerized modula- 
tion of hopping amplitudes as depicted in Fig. 1(a). Such 
modulations can arise as a Peierls distortion of a uniform ion 
lattice above a critical value of the electron-phonon coupling 
or, as discussed in Ref. |4| in the context of graphene, as an 
interaction-driven instability. We show that a vortex in the 
complex scalar order parameter describing this dimerization 
pattern generates a zero-energy bound state in the spectrum of 
electrons and carries a fractional charge. In addition to simi- 
lar arguments to those of Ref. |4j, we present a simple electron 
counting argument 5,6 as well as numerical evidence demon- 
strating this effect. Importantly, we find that the fractionaliza- 
tion survives essentially intact beyond the low-energy theory 
of Ref. 4 in the presence of the lattice. We also show that the 
same pattern of fractionalization occurs for the Z4 and U(l) 
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FIG. 1: (Color online) The model: (a) Square lattice with ^Oo mag- 
netic flux per plaquette and dimerized hopping amplitudes. The + on 
the left for each row show the choice of gauge for the Peierls phase 
factors. The solid (dashed) bonds indicate an increased (decreased) 
hopping amplitude in the x (blue) and y (green) directions as ex- 
plained in the text. The four sites of the unit cell are marked, (b) The 
Z4 vortex. The dashed lines indicate the domain walls sharing the 
center of the vortex. The phase of the local dimer order parameter f 
(see text) around (c) the U(l) and (d) the Z 4 vortex. The x shows the 
center of each vortex where / = 0. 



vortices and clarify the issue of their confinement. We then 
give a brief discussion of the energetics of the Peierls distor- 
tion and touch upon possible experimental realizations of the 
model in artificially engineered semiconductor heterostruc- 
tures and optical lattices. 

The model. — We consider a square lattice with a tight- 
binding Hamiltonian 
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for electrons annihilated at site r, by c, . We assume the spins 
of all electrons are polarized along the field. Here, fy are hop- 
ping amplitudes between nearest-neighbor sites and e, is the 
on-site potential. The magnetic field is included through the 
Peierls phase factors, 0y = (2^/<D>o) j£ 1 A ■ d\. Each plaquette 
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is threaded by flux <I>o/2. We choose to work in the Landau 
gauge A = (<l>()/2)(-i/, 0) where we have set the lattice spac- 
ing to unity. The hopping amplitudes in our model form a 
dimerized pattern as shown in Fig.[Tf a) with alternating bonds 
tij+ir = f(l ± m x ) and ?,•_,-+§ = f(l + m y ). In the following we 
assume uniform e; = e. 

We can now arrange the four sites of the unit cell into a 
spinor field 

*Ak = e^ 72 (ci k , c 2k , c 3k , c 4k ) r (2) 
to write the Hamiltonian as 

We = J] <A k (e-2f#k)<Ak (3) 

keBZ 

in the reduced Brillouin zone BZ = {k||JU,|iy < \n} with 
#k = To (n cos + y 2 cos k y + m x sin K + im y y 5 sin fe^J . 

We are using the standard Weyl representation of Dirac matri- 
ces: yo — o~\ <8> 1, y — icr 2 ® c? and 75 = -iyo7i7273 = 0"3 ® 1. 
The spectrum of H k is given by 



£ k = ± -ycos 2 fc x + cos 2 k y + m\ sin 2 k x + m 2 sin 2 (4) 

which is doubly degenerate and symmetric: if if/ E is an eigen- 
state with energy E, i.e. H^e = Eif/£, then yoy^E is an 
eigenstate with energy —E since {H k , 7073} = 0. 

Low-energy theory. — The undistorted system with m = 
(m x ,m y ) = has Dirac nodes at k x — k y — +n/2 where 
the two energy bands meet. There is one node per reduced 
Brillouin zone. At nonzero constant dimerization m a gap 
m = (m 2 + ffJ 2 ) 1 ^ 2 opens up in the spectrum. We will focus our 
attention on the low-energy excitations of the system in the 
situation when the dimerization vector m is slowly varying in 
space. Using a polar representation for the hopping anisotropy 
m x + irriy = me lx and by linearizing 'He around the nodes in 
the continuum limit we find the low-energy Hamiltonian 

H x = 70 (in d x + iy 2 d y + me lxj5 ) . (5) 

An interesting situation arises when a U(l) vortex is present 
in the dimerization vector m. Such a vortex is characterized 
by a 2nn winding in the phase x> with n an integer, along 
a contour enclosing the vortex center. The continuum Dirac 
Hamiltonian <(5j in this vortex background field is the same as 
that studied in Refs. Iand0 It admits \n\ degenerate zero- 
energy bound states. 8 This result has its origin in an index 
theorem for the Dirac Hamiltonian in topologically nontrivial 
background fields. 9 

For a singly quantized vortex (|ra| = 1) there is a single zero 
mode with support on one of the sublattices depending on the 
sense of the vortex. A standard calculation^ - can be used to 
show that when the lattice is half-filled the charge carried by 
such a vortex is fractionalized, namely ±e/2, depending on 
whether or not the zero mode is occupied by an electron. 

In our square lattice model one can also give a simple intu- 
itive counting argument for the charge fractionalization that 



generalizes the construction deployed in the case of poly- 
acetylene by Goldstone and Wilczek. 5 Consider a Z4 vortex, 
depicted in Fig. l(b,d), in the limit of strong dimerization, 
m =s 1. In this limit each thick bond in Fig. 1(b) can be 
thought of as representing one electron resonating between 
the two lattice sites. The topological structure of this Z4 vor- 
tex guarantees that in its core there is necessarily a site with 
no partner. 1 1 This site can be either occupied or empty, thus 
exhibiting e/2 surplus or deficit of charge compared to the 
background, which has e/2 per site on average. The Z4 vor- 
tex can be adiabatically deformed into the U(l) vortex of Fig. 
1(b) and to the weakly dimerized phase at m <k 1. This indi- 
cates that the charge fractionalization persists in that limit, as 
explicitly verified by the analytical calculation 410 and by the 
numerical results reported below. We note that this argument 
works in our model only in the presence of the magnetic flux, 
which is needed for the existence of the Dirac nodes in the 
spectrum. At zero flux there is a Fermi surface for any small 
m and the unpaired center of a vortex does not create a bound 
state. 

Numerics. — We have performed exact diagonalization of 
the tight-binding Hamiltonian ([TJ with different vortex con- 
figurations on lattices of size up to 52 x 52. Open boundary 
conditions were assumed in all cases. We consider two re- 
alizations of the vortex, a Z4 vortex and a lattice version of 
the continuum U(l) vortex studied in the low-energy theory. 
To quantify the vortex structure we use a local dimer order 
parameter 12 defined at site r, = (x,, yi) as 

fi = ^ nu+fiki+P- ( 6 ) 

Here the fi are the four nearest-neighbor unit vectors, T]ij +X = 
(-1)*', T] U+ g = i(-l) y > and 77,^ = tj^. (See also Fig. 2 in 
Ref.0) 

The Z4 vortex occurs if the energetics of the problem allow 
only dimerizations along the high-symmetry directions of the 
crystal, e.g. m = (+m, 0) and (0, +;«). A vortex is then ob- 
tained by putting together four possible domains as illustrated 
in Fig. \V[b), separated by domain walls emanating from the 
vortex center. The phase of the order parameter for the Z4 
vortex is depicted in Fig. |TJd). 

The U( 1 ) vortex occurs if the energy of the state character- 
ized by dimerization vector of the form m = m(cos^-, sin^-) is 
approximately independent of x- We construct the U(l) vor- 
tex on the lattice by discretizing the continuum vortex field 
(p(z) = z/\z\ where z — x + iy is the complex coordinate of the 
point (x, y). Explicitly, 

fi,;+p = fo(l + mRe [vti +fi <P (t { + jfi)] ), (7) 

where the asterisk denotes the complex conjugation. Fig.QJb) 
shows the phase of /; clearly illustrating the U(l) form. The 
U(l) structure can be extended to study the case for multiple 
vortices and antivortices. For vortices atzi, Zz, ■ ■ ■ and antivor- 
tices at z[ , z' 2 , ■ ■ ■ this is achieved by taking a product of vortex 
fields <p(z) = Ua(z - z a )/\z - z a \ Ub(z - z' b Y l\z - z' b \ in Eq. 0. 

A typical sample of our results for a system with one and 
two vortices is shown in Figs.|2]and[3] respectively. The struc- 
ture of the vortex is reflected in the symmetries of the zero 
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FIG. 2: (Color online) The charge density at half-filling with zero 
mode filled in the U(l) vortex with m = 0.15 (a,b). The same plots 
for the Z4 vortex and m = 0.30 (c,d). The lattice size 51x51 and 
the gridline spacing is 3. The top- view density plots on the left are 
included to show the symmetry of the solution for each vortex and 
the support on a single sublattice. 
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FIG. 3: The charge density at half-filling with both midgap states 
filled (a,c) and the low-energy (including midgap) eigenvalues (b,d) 
of a system with two U(l) vortices (a,b) and with a vortex and an 
antivortex (c,d). The gap m = 0.2, the lattice size is 52 x 51 and the 
gridline spacing is 3. 



mode and the corresponding charge density at half-filling. For 
a single vortex, Fig. [2] the charge density at half-filling with 
the zero mode filled shows an excess charge localized at the 
vortex above the background value. For both the U(l) and Z4 
vortex this excess charge integrates to 0.5e to within machine 
accuracy confirming the existence of the fractional charge. 
Remarkably, the charge density of the half-filled system with 
the Z4 vortex is essentially the same as that with the U(l) vor- 
tex. Especially, no charge is bound to the domain walls far 
away from the center. 

The system with two vortices has two zero-energy states 
with support on the same sublattice, Fig. [3ja). It is impor- 
tant to note that, as illustrated in Fig.[3jb), these remain exact 
zero modes even when the intervortex separation d vv becomes 
short compared to the size of the bound state wavefunction 
£ ~ m~ x . The zero modes are topologically protected in accor- 
dance with the index theorem of Ref. 9 which states that the 
number of exact zero modes in the system is bounded from 
below by the total vorticity of the order parameter /. This can 
be understood on a more intuitive level by noting that since 
the two bound states reside on the same sublattice the matrix 
element of the Hamiltonian ([TJ with e — between the two 
states vanishes identically and hence leads to no splitting in 
energies. Alternatively, one can consider a process of merging 
n single vortices into one vortex with winding n. According 
to the Ref.[||the latter supports \n\ zero modes. Thus, the zero 
modes associated with n individual widely separated vortices 
transform smoothly into n zero modes of the n-fold vortex. 

The vortex-antivortex configuration, by contrast, does not 
have exact zero modes. Now the two bound states reside 
on different sublattices, Fig. |3jc), so there is a finite tunnel- 
ing amplitude between them. This yields two midgap states 
with energies +E m that depend exponentially on the vortex- 



antivortex separation d vd and merge with the continuum states 
as d va — > 0, Fig. [3jc). Again, these results are consistent 
with the index theorem 9 since the total vorticity of a vortex- 
antivortex pair vanishes. 

Peierls distortion. — In a crystalline system the dimeriza- 
tions postulated in our model may occur through a Peierls dis- 
tortion of the ion lattice in the presence of the magnetic flux at 
or near half-filling. Such a distortion is driven by the electron 
kinetic energy gain upon the opening of a gap. We briefly out- 
line this mechanism. At half filling, the energy gain per site is 
5& e = — — 2k[£k(m) - £k(0)], where N is the number of sites 
and £k(ni) is given by the positive branch of Eq. ©. This gain 
is to be balanced against the loss of elastic energy that results 
from the ion displacements u,, 6&\ m = jb 2<i/>( u i _ u j) 2 - We 
assume henceforth that the change in the hopping element is 
related to the displacement as m l ;+ p = gfi ■ (u, - u, + /j) with g 
the dimensionless electron-phonon coupling. 13 For the dimer- 
ization pattern characterized by a constant, spatially uniform 
vector m, we have d>£i att /4f = i/em 2 with k = K/4tg 2 . 

The Peierls distortion occurs when 6S = 6S e + ££i att be- 
comes negative. Unlike the one-dimensional case, where the 
Peierls instability occurs at infinitesimal g, in two dimensions 
g must exceed a critical value. In terms of k this happens for 
K < k c = N~ l 2k sin 2 k x /E k (0) - 0.806. Near the transition 
we can expand 

-77 = » (* - KcW + — M 3 + 0(m 4 ). (8) 

4f 2 J7T 

We observe that the leading terms in 6S are isotropic in m. 
The anisotropy first enters at the fourth order via a term 
~ (m\ - m 2 ) 2 . This indicates that for small m a vortex in this 
model will closely resemble the U(l) vortex in Fig. 1(c), with 
a small four-fold anisotropy. For m 1 the anisotropy be- 
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comes important and the pattern will be closer to a Z4 vortex, 
Fig. 1(d). 

These findings are important since the Z4 vortices are gener- 
ically confined by a linear potential due to the domain walls 
while the U(l) vortices are not. For the fractional charges to 
be observable we require that the Z4 confinement length scale 
be much larger than It is also possible that the U(l) sym- 
metry emerges universally at the dimerization transition. We 
note that in this case the small Z4 anisotropy of the present 
model is an irrelevant perturbation while the Z3 in graphene 4 
is relevant. 

Experimental realization. — In a natural solid with lattice 
constant arj ~ 2A the magnetic field required to achieve the 
Oo/2 flux per plaquette is of the order 10 4 T and thus much 
too high to reach in a laboratory. A better chance for observ- 
ing the phenomena discussed here is to artificially engineer 
the hosting system such that the electrons move in a lattice 
with much larger ao- 

Indeed several schemes have been proposed and imple- 
mented recently to achieve this effect in semiconductor het- 
erostructures (such as GaAs/AlGaAs) combined with a type-II 
superconducting layer 3,14 or a top gate realized by the diblock 
copolymer nanolithography technique^ A second possibility 
is offered by fermionic cold atoms in an optical lattice with 
artificial magnetic flux. 16 Once the tight binding model with 
the magnetic flux is realized by any of the above techniques 
there of course remains the challenge of creating a vortex and 
detecting the fractional charge. 



Conclusion. — We showed through general arguments and a 
detailed numerical study that the system proposed here sup- 
ports configurations with fractional charge while preserving 
time-reversal symmetry. The number of zero modes in vortex 
configurations leading to this fractionalization is topologically 
protected. This should be contrasted with similar zero modes 
in vortices of a chiral p-wave superconductor where such a 
generic topological protection seems not to exist^ An in- 
teresting question is what exchange statistics these fractional 
particles obey. This is significant theoretically as well as for 
potential applications in quantum information processing. A 
natural expectation, based on the analogy with the fractional 
quantum Hall states, 1,2 is that they are anyons. This is also 
suggested by arguments using gauge invariance on a torus. 18 
Since the system is time-reversal invariant one expects the ef- 
fective theory for the topological excitations to be a "doubled" 
Chern-Simons gauge theory of the type discussed in Ref. [l9[ 
Indeed, our preliminary investigation 20 suggests that such a 
gauge structure will emerge upon integrating out the Fermi 
fields in the Dirac Hamiltonian (0 in the presence of topolog- 
ically nontrivial background field \- 
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